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Abstract. - We use a density- functional-theory (DFT) approach with a modified Becke-Johnson 
exchange plus local density approximation (LDA) correlation potential (mBJLDA) [semi-local, 
orbital-independent, producing accurate semiconductor gaps, see F. Tran and P. Blaha, Phys. 
Rev. Lett. 102, 226401 (2009)] to investigate the electronic structures of zincblende transition- 
metal (TM) pnictides and chalcogenides akin to semiconductors. Our results show that this 
potential does not yield visible changes in wide TM d-t2 g bands near the Fermi level, but makes 
the occupied minority-spin p-bands lower by 0.25~0.35 eV and the empty (or nearly empty) 
minority-spin e g bands across the Fermi level higher by 0.33^0.73 eV. Consequently, mBJLDA, 
having no atom-dependent parameters, makes zincblende MnAs become a truly half-metallic (HM) 
ferromagnet with a HM gap (the key parameter) 0.318eV, being consistent with experiment. For 
zincblende MnSb, CrAs, CrSb, CrSe, or CrTe, the HM gap is enhanced by 19~56% compared to 
LDA and generalized gradient approximation results. The improved HM ferromagnetism can be 
understood in terms of the mBJLDA-enhanced spin exchange splitting. 
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Introduction. — It is believed that half-metallic 
ferromagnetic materials can be used to make high- 
performance spintronic devices [T][2] , because one of their 
two spin channels is metallic and the other has the fea- 
ture of a semiconductor or insulator, and then the spin 
polarization at the Fermi energy can reach 100% if the 
spin-orbit coupling effect is neglected [3J. Half-metallic 
ferromagnetism has been found in NiMnSb [3J, Cr02 @], 
?e 3 4 0, Co 2 MnSi 0, and La .67Sro.33Mn0 3 0, and 
nearly 100% high spin-polarization has been observed ex- 
perimentally in the cases of CrC>2 and Lao.67Sro.33Mn03 
materials [8]. It is of much interest to seek half-metallic 
ferromagnetic materials compatible with important semi- 
conductors. Great efforts have been made in this di- 
rection, and significant progresses have been achieved 
in binary transition-metal pnictides and chalcogenides 
such as MnAs, MnSb, CrAs, CrSb, CrSe, and CrTe 
with zincblende [9TI25] and wurtzite structures [26], spe- 
cial transition- metal-doped semiconductors [23HH]: other 
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transition-metal compounds based on convectional semi- 
conductors |2"9T[3Tj etc. 

On the theoretical side, one usually uses density- 
functional-thcory (DFT) [32] method with generalized 
gradient approximation (GGA) [33], local density approx- 
imation (LDA) [31] etc to study the electronic structures 
of these materials akin to semiconductors. Because it is 
well-known that GGA and LDA underestimate semicon- 
ductor gaps, there naturally exists a question: How good 
do GGA and LDA describe the electronic structures of 
such half-metallic ferromagnetic materials? This needs 
to be answered because there exists some inconsistence 
between theoretical predictions and corresponding experi- 
mental observations. For examples, zincblende MnAs has 
been experimentally shown to be half- metallic [TQH]~2] , but 
both GGA and LDA calculations prove that it is not truly 
half-metallic [3] because the Fermi level touches or even 
crosses the bottom of the minority-spin conduction bands. 
Fortunately, there have already been many feasible ap- 
proaches to improve DFT calculations of semiconductor 
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gaps, such as various GW methods [551156] . exact-exchange 
approach [3"7H3"9"] . different hybrid- functionals [ID], and a 
modified Becke- Johnson (mBJ) exchange potential [3T] . 
Some of them demand high expenses in computational 
resources, and some are very economical. While be- 
ing as computationally economical as GGA and LDA 
and needing no atom-dependent parameters, the mBJ ex- 
change plus LDA correlation (mBJLDA) potential as an 
orbital- independent, semi- local exchange-correlation po- 
tential has been proved to produce accurate gaps for wide- 
band-gap insulators, sp semiconductors, and 3d transition- 
metal oxides [41] . Therefore, it is highly desirable, even 
necessary, to use this greatly-improved, atom-parameters- 
independent, still economical DFT approach to investigate 
such half-metallic materials akin to semiconductors or in- 
sulators. 

In this article we use a state-of-the-art DFT method 
with mBJLDA potential [3JJ to investigate the electronic 
structures of half-metallic ferromagnets akin to semicon- 
ductors. Our results show that in the cases of zincblende 
MnAs, MnSb, CrAs, CrSb, CrSc, and CrTe, mBJLDA 
does not yield visible changes in wide transition-metal 
d-^2g dominating bands near the Fermi level, but makes 
the occupied minority-spin p-dominating bands lower by 
0.25—0.35 eV and the empty (or nearly empty) minority- 
spin e g bands higher by 0.33—0.73 cV compared to 
GGA and LDA results. Consequently, mBJLDA makes 
zincblende MnAs become a truly half-metallic ferromag- 
net, being consistent with experiment, and enhances the 
half-metallic gaps [32] of the other five by 19-56%. The 
improved half-metallic ferromagnetism can be understood 
in terms of the fact that the spin exchange splitting is 
enhanced by using mBJLDA potential. This consistence 
is a direct evidence that mBJLDA as a semi-local and 
orbital-independent exchange-correlation potential with- 
out atom-dependent parameters is very satisfactory in de- 
scribing the electronic structures of such materials akin to 
semiconductors. More detailed results will be presented 
in the following. 

The rest of the paper is organized as follows. We shall 
describe our computational detail in next section. In the 
third section we shall present our main results calculated 
with mBJLDA in comparison with those calculated with 
the GGA and LDA. In the fourth section we shall make 
necessary discussions and try to understand the calculated 
results. Finally, we shall give our conclusion in the fifth 
section. 

Computational detail. — We use a full-potential lin- 
earized augmentcd-plane-wavcs method within the den- 
sity functional theory [32] . as implemented in package 
WIEN2k [43] . We use mBJLDA for the exchange- 
correlation potential [34l[4T] to do our main DFT calcula- 
tions, and take popular LDA by Perdew and Wang (LDA- 
PW91) [3JJ and GGA by Perdew, Burke, and Ernzerhof 
(GGA-PBE) [231 to do corresponding comparative stud- 
ies. The mBJ exchange potential [3T] was developed from 



a semi-local exchange potential proposed by Becke and 
Johnson (BJ) [44] . The BJ exchange potential can repro- 
duce the shape of the exact-exchange optimized-effective 
potential of atoms [44]. mBJLDA as a combination of 
mBJ exchange and LDA correlation can produce accurate 
semiconductor gaps for sp semiconductors, widc-band-gap 
semiconductors, and transition-metal oxide semiconduc- 
tors and insulators [31]. The full rclativistic effects are 
calculated with the Dirac equations for core states, and 
the scalar relativistic approximation is used for valence 
states [35H37] . The spin-orbit coupling is neglected be- 
cause it has little effect on our results. We use 3000 k 
points in the first Brillouin zone, make harmonic expan- 
sion up to Z max = 10 in each of the atomic spheres, and 
set R m t * fcmax = 8. The radii of the atomic spheres of 
Mn, Cr, As, Sb, Sc, and Tc are set to 2.41, 2.41, 2.10, 
2.27, 2.16, and 2.33 bohr, respectively. The volumes are 
optimized by using LDA-PW91 and GGA-PBE in terms 
of total energy method. The self-consistent calculations 
are considered to be converged when the integration of 
absolute charge-density difference between the input and 
output electron density is less than 0.0001 |e| per formula 
unit, where e is the electron charge. 

Main calculated results. — At first, we investigate 
zincblende MnAs to show advantage of mBJLDA. This is 
because on the experimental side it has been shown to 
be half-metallic and intensively studied for device appli- 
cations [T0HT2] , but on the theoretical side both LDA and 
GGA calculations show that it is not truly half- metallic [5] . 
Since there is no corresponding energy functional for the 
mBJ exchange potential [41j . we use experimental lattice 
constant as our computational input. We study its elec- 
tronic structures by using the three exchange-correlation 
potentials: LDA-PW91, GGA-PBE, and mBJLDA. We 
present its total and partial DOSs calculated with mB- 
JLDA, GGA-PBE, and LDA-PW91 in Fig. 1. The mBJ 
result shows that zincblende MnAs is a half-metallic ferro- 
magnet. The total magnetic moment per formula unit is 
4/iB, in contrast with 3.86/is calculated with GGA-PBE. 
The minority-spin gap across the Fermi level is 2.352 eV, 
substantially larger than the GGA-PBE result, 1.76 cV, 
and the half-metallic gap is determined by the bottom 
of the minority-spin conduction bands, being 0.318 eV. 
Similar results could be obtained by adding appropriate 
orbital-dependent on-site correlation parameters to some 
localized orbitals, as was done in LDA+U approach [48] 
in the case of Co2FeSi [4jJ[. Another special feature in 
the mBJLDA DOS is that the As-p-dominating bands are 
merged with the Mn-d-dominating bands in the majority- 
spin channel, and the As-p-dominating bands in both spin 
channels are a little narrower than those calculated with 
both GGA-PBE and LDA-PW91. The mBJLDA DOS 
result implies 100% spin polarization at the Fermi level, 
but either LDA-PW91 or GGA-PBE result leads to much 
smaller spin polarization at the Fermi level. Therefore, the 
mBJLDA DOS is consistent with experiment, but either 
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of the LDA-PW91 and GGA-PBE results is not. 
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Fig. 1: (color online) Spin-resolved density of states of MnAs 
with the zincblende structure, calculated with mBJLDA (up- 
per panel), GGA-PBE (middle panel), and LDA-PW91 (lower 
panel), respectively. The upper part in each panel is for 
majority-spin channel and the lower part for minority-spin. 



In Fig. 2 we present the spin-resolved energy bands be- 
tween -5 and 4.5 eV of zincblende MnAs calculated with 
mBJLDA and GGA-PBE. There are little differences be- 
tween mBJLDA and GGA for the majority-spin t2 g bands 
near the Fermi level, but the mBJLFA raises the minority- 
spin e g bands by 0.3 eV and lowers the majority-spin e g 




Fig. 2: (color online) Spin-resolved energy bands of MnAs with 
the zincblende structure, calculated with mBJLDA (solid lines) 
and GGA-PBE (dashed lines), respectively. The left panel is 
for majority-spin channel and the right panel for minority-spin. 




Fig. 3: (color online) The charge (upper panel) and magne- 
tization (lower panel) density distributions of MnAs with the 
zincblende structure. The plane is defined by an As-Mn-As 
triangle. Upper: The actual charge density is 0.004~0.068 
|e|/bohr 3 and the density increment is 0.004|e|/bohr 3 . Lower: 
The lines with '0' label the zero value for the magnetiza- 
tion density (-0.008^+0.008 |^s|/bohr 3 ) and the increment 
is 0.0008/ig/bohr 3 . 



bands by 0.3 eV. These are vital to make zincblende MnAs 
become a truly half-metallic ferromagnet. We present 
the charge and magnetization density distributions of the 
zincblende MnAs calculated with the mBJLDA potential, 
within the charge density range - 0.068 |e|/bohr 3 and 
the magnetization density -0.008 ~ +0.008 l^el/bohr 3 , 
respectively. Here the energy range is from -5 eV to the 
Fermi level and the plane is defined by the nearest As-Mn- 
As triangle. There is not any visible difference if mBJLDA 
is replaced by GGA. It can be seen that the Mn-As bond 
is covalent and there is a small negative magnetization 
around As although the total magnetic moment in the 
cell is contributed mainly by Mn. 

The mBJLDA potential has similar effect on the elec- 
tronic structures of other zincblende compounds, such as 
MnSb, CrAs, CrSb, CrSe, and CrTe. We present the spin- 
resolved total DOSs of zincblende MnAs, MnSb, CrAs, 
CrSb, CrSe, and CrTe calculated with GGA-PBE and mB- 
JLDA. For convenience of comparison, we plot the GGA- 
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Fig. 4: (color online) Spin-resolved density of states of MnAs (left-upper), MnSb (left-lower), CrAs (middle- upper), CrSb 
(middle- lower), CrSe (right-upper), and CrTe (right-lower) with the zincblende structure, calculated with mBJLDA (solid lines) 
and GGA-PBE (dashed lines), respectively. The upper part in each panel is for majority-spin channel and the lower part for 
minority-spin. 



PBE and mBJLDA DOSs in the same panel for each of 
the six zincblende compounds. It is clear that in each of 
the six cases the minority-spin gap across the Fermi level 
is made larger and the majority-spin e g bands move down- 
ward when GGA-PBE is replaced by mBJLDA. Further- 
more, the bottom of the minority-spin conduction bands 
calculated with mBJLDA is higher than that with GGA, 
and consequently, the half-metallic gap is effectively made 
larger by replacing GGA with mBJLDA. All the mBJLDA 
results of the magnetic moments (M), the minority-spin 
gaps (Gmis)j and the half-metallic gaps (Ghm) are sum- 
marized in Table [TJ For comparison, we also present 
the corresponding GGA-PBE results in the parentheses. 
It is clear that both G m i s and Ghm are made larger by 
replacing GGA-PBE with mBJLDA. Most importantly, 
zincblende MnAs is a truly half-metallic fcrromagnet if 
mBJLDA is used. Therefore, the half-metallic fcrromag- 
nctism has been improved after we use mBJLDA instead 
of GGA-PBE. 

Discussions. — It has been proved that mBJLDA 
potential improves semiconductor gaps for various semi- 
conducotrs [U]. Our calculated results have shown that 
the minority-spin gap across the Fermi level and the 
half- metallic gap (the key parameter [42]) of the binary 
transition-metal pnictides and chalcogenides are both im- 
proved by mBJLDA potential. Especially, the nearly- 
half-metallic zincblende MnAs under GGA and LDA be- 
comes a truly half-metallic ferromaget when we use mB- 
JLDA potential. This is consistent with experimental re- 



Table f : The lattice constants (a), magnetic moments per for- 
mula unit (M), minority-spin gaps (Gmis), and half-metallic 
gaps (G hm ) of zincblende MnAs, MnSb, CrAs, CrSb, CrSe, 
and CrTe calculated with mBJLDA. The values in the paren- 
theses are calculated with GGA-PBE. 



Name 


a (A) 


M {fx B ) 


G mis (eV) 


G hm (eV) 


MnAs 


5.717 


4 (3.8) 


2.352 (1.760) 


0.318 ( - ) 


MnSb 


6.166 


4(4) 


1.923 (1.491) 


0.459 (0.200) 


CrAs 


5.659 


3(3) 


2.579 (1.896) 


0.636 (0.322) 


CrSb 


6.138 


3(3) 


2.180 (1.650) 


0.992 (0.751) 


CrSe 


5.833 


4(4) 


4.569 (3.508) 


1.331 (0.604) 


CrTe 


6.292 


4(4) 


3.692 (2.919) 


1.384 (0.858) 



suit that zincblende MnAs is a truly half-metallic ferro- 
magnet [TTJHT2"] . Therefore, we believe that mBJLDA po- 
tential as an exchange-correlation potential is more sat- 
isfactory than both GGA and LDA in describing half- 
metallic ferromagents akin to corresponding semiconduc- 
tors, such as zincblende MnAs, MnSb, CrAs, CrSb, CrSe, 
and CrTe. More fortunately, the computational efficiency 
of mBJLDA is almost the same as those of GGA and LDA, 
and hence one can calculate large systems with mBJLDA 
as one does with GGA and LDA. 

What causes this improvement? On one hand, our cal- 
culated results show that the replacing GGA with mB- 
JLDA has little effect on the wide majority-spin bands 
near the Fermi level. Actually, mBJLDA does not yield 
visible changes in wide transition-metal d-t2 g bands with 
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small density of states, but indeed shifts the narrow e g 
bands with respect to the t^g bands. We can see in Fig. 
4 that the occupied majority-spin e g bands emerge with 
the lower occupied majority-spin p-bands in the cases 
of MnAs and MnSb. On the other hand, the occupied 
minority-spin p-bands arc lowered by 0.25—0.35 cV and 
the empty (or nearly empty) minority-spin e g bands are 
raised by 0.33—0.73 eV, which enhances substantially both 
the majority-spin gap and the half-metallic gap. These 
can be understood by assuming it is the mBJLDA po- 
tential that enhances the spin exchange splitting with 
respect to GGA and LDA results. This explanation is 
supported by the observation in Fig. 1 that there are 
much d-t2 g weight in the top of the minority-spin valence 
bands and in the bottom of the minority-spin conduction 
bands. The enhancement of the spin exchange splitting 
is achieved because of the advantage of the BJ exchange 
potential [5TJ[S]. Thus, the half- metallic ferromagnctism 
is enhanced in the binary transition-metal pnictides and 
chalcogcnides akin to semiconductors. 

Conclusion. — In summary, we have used a state- 
of-the-art DFT method with the mBJLDA potential [MJ 
|4"T] to investigate the electronic structures of zincblende 
transition-metal pnictides and chalcogcnides and compare 
the mBJLDA results with those calculated with popu- 
lar exchange-correlation potentials (GGA-PBE and LDA- 
PW91). Our calculated results show that the mBJLDA 
potential does not yield visible changes in wide transition- 
metal d-t2 g bands near the Fermi level, but makes the oc- 
cupied minority-spin p-bands lower by 0.25—0.35 eV and 
the empty (or nearly empty) minority-spin e g bands higher 
by 0.33—0.73 eV. In the cases of zincblende MnSb, CrAs, 
CrSb, CrSe, and CrTe, the half-metallic gaps are enhanced 
by 19-56% with respect to the GGA-PBE and LDA- 
PW91 results. Most importantly, the mBJLDA poten- 
tial makes zincblende MnAs become a truly half-metallic 
ferromagnet. This result is consistent with experiment, 
which is a direct evidence that the mBJLDA potential 
is more satisfactory than popular GGA and LDA in de- 
scribing the electronic structures of the binary transition- 
metal pnictides and chalcogcnides akin to semiconductors. 
The improved half-metallic ferromagnctism can be natu- 
rally understood in terms of the enhanced spin exchange 
splitting caused by the mBJLDA potential. We believe 
that the semi-local and orbital-independent mBJLDA po- 
tential without atom-dependent parameters will be more 
satisfactory than usual GGA and LDA in calculating elec- 
tronic structures with almost the computational efficiency 
of usual GGA and LDA for magnetic materials sharing 
some features with semiconductors. 
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